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Abstract. In 1973, Lc Bellac and Levy-Leblond (Nuovo Cimento B 14 217-234) 
discovered that Maxwell's equations possess two non-relativistic Galilei-covariant 
limits, corresponding to |E| ^> c|B| (electric limit) or |E| <C c|B| (magnetic limit). 
Here, we provide a systematic, yet simple, derivation of these two limits based on 
a dimensionless form of Maxwell's equations and an expansion of the electric and 
magnetic fields in a power series of some small parameters. Using this procedure, 
all previously known results are recovered in a natural and unambiguous way. Some 
further extensions are also proposed. 
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1. Introduction 

Maxwell's equations were the first physical theory to be put forward that is fully Lorentz 
covariant, well before the special theory of relativity was developed and understood. 
Indeed, it was the clash between the Lorentz covariance of Maxwell's equations and 
the Galilei covariance of Newtonian mechanics that stimulated the discovery of special 
relativity and revealed that Lorentz covariance is the most fundamental symmetry of 
the two. 

Nevertheless, it should be stressed that electricity and magnetism were developed 
historically as independent phenomena and only lately were realized to be different 
expressions of a single underlying theory. Already in the eighteenth century, physicists 
like Charles-Augustin Coulomb (1736-1806) manipulated electric charges and measured 
how they attract or repel each other through electric fields - in modern language, they 
worked out applications of Gauss's law: V • E = p/eq. This is the field of electrostatic. 

One century later, physicists working on magnetism manipulated currents and 
measured how they interact with each other through magnetic fields. Around 1820, 
Oersted observed that wires carrying electric currents deflected a compass needle placed 
in their vicinity. Biot and Savart, and later Ampere, established rigorous laws that 
related the strength and direction of a magnetic field to the currents that produce it. In 
modern language, they worked out applications of Ampere's law: V x B = /ioJ. This 
is the field of magnetostatics. 

But none of these phenomena involve anything that is "relativistic" . This is obvious 
for electrostatics, as only electric charges and fields are involved; but it is true for 
magnetostatics too, because the current in Ampere's law does not need to be relativistic 
in order to generate a finite magnetic field. 

Thus, we are faced with two perfectly valid theories that appear to be both 
non-relativistic and both deriving from the complete theory of electromagnetism, i.e. 
Maxwell's equations. But surely any relativistic theory should have a unique non- 
relativistic limit - or should it? 

Enter Le Bellac and Levy-Leblond (hereafter, LBLL), who in 1973 published a 
paper where it was demonstrated that Maxwell's equations possess two distinct, yet 
perfectly consistent, non-relativistic limits pp. The two limits correspond to situations 
where either |E| ^> c|B| (electric limit) or |E| <C c|B| (magnetic limit). Each of the 
two limits is Galilei covariant, although the transformations of the fields and of the 
density and current are not the same in the two cases. In practice, the electric limit 
amounts to neglecting the time-derivative of the magnetic field in Faraday's law of 
induction, whereas the magnetic limit is obtained by dropping the displacement current 
in Ampere's equation. 

LBLL's paper basically contains all that one needs to know on this topic, but also 
has several drawbacks. First, the results were presented without a systematic derivation; 
only a posteriori it was checked that both limits are indeed Galilei covariant. Second, the 
authors obtained their result by employing SI units, which was sort of a novelty at their 
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time of writing, CGS Gaussian units being then the preferred choice. But any physically 
meaningful result should not depend on the units in which the equations are written, 
and a proper limit should be obtained by making use of dimensionless quantities. 

More recently, several papers revisited the work of LBLL from different angles, 
elucidating some issues such as gauge invariance [2] and the correct Galilean limit of the 
relativistic four-current [3], discussing applications [I], or else extending the analysis to 
the coupled Dirac-Maxwell equations [5j. Some papers use more abstract methods that 
rely on a five- dimensional Minkowski manifold [HJ [7]. 

In none of these recent works the electric and magnetic limits were derived 
systematically from the full Maxwell equations. In general, the relations |E| 3> c|B| 
or |E| c|B| (and the analogue expressions for the sources) are assumed ex nihilo and 
are used to derive the low-velocity transformations of the fields and sources. It is then 
proven that some form of "reduced" Maxwell equations are Galilei invariant under those 
transformations [5j. Only in one case this analysis is carried out using dimensionless 
variables [8] (see also [9]). Sometimes, non-systematic "order of magnitude" arguments 
are used heuristically to justify the limits [TO] . 

The purpose of the present paper is to introduce a simple, yet systematic, procedure 
to recover the non-relativistic limits of Maxwell's equation. This procedure is based 
on Maxwell's equations written in a non-dimensional form. It is shown that two 
independent dimensionless parameters naturally appear in the equations. The electric 
and magnetic limits are then derived by letting either or both these parameters go 
to zero. Using the same procedure, all known properties (gauge relations, Lorentz 
transformations, the Lorentz force) of the non-relativistic limits can be deduced 
systematically. This procedure is also better adapted to possible extensions of the 
present work, such as the coupling of the Maxwell and Dirac equations in relativistic 
mean field theories [TTl [T2] . 



2. General procedure 

We start from Maxwell's equations in SI units: 

VE =A (1) 

V-B =0, (2) 
<9B 

VxE=--, (3) 
1 <9E 

VxB = ^ J + --, (4) 

& at 

where Eq and fiQ are respectively the electric permittivity and the magnetic permeability 
in vacuum, and eo/^o = c ~ 2 - 

We normalize space to a reference length L and time to a reference time T, which 
define a typical velocity V = L/T. The fields and the sources are also normalized to 
reference quantities denoted by an overline: E, B, ~p, and J, satisfying the relations: 
E = cB and 7 = Vp. 
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In these units, Maxwell's equations can be rewritten as: 

V-E (5) 

a 

VB =0 (6) 

VxE=-*f (7) 

VxB = £ J + /J f, (8) 
where the following two dimensionless parameters have appeared naturally: 

/? = -, a=|Si. (9) 
c pL 

The first of these parameters, /3, is just the reference velocity normalized to the speed 
of light and is obviously small in any non-relativistic limit. 

The parameter a requires a little analysis to be expressed in more familiar terms. 
We can think of our distribution of charges and currents as a classical plasma, with 
a certain temperature T and typical charge density p. We then define a reference 
electrostatic potential = EL and express it in terms of the temperature, such that 
qcf) = ksTo, where q is the electric charge of the particles and ks is the Boltzmann 
constant. Then a can be written as: 



a = ( -f ) (10) 



where Xo = \Zk~BToSoJ{qpj is the Debye length. In plasma physics, this ratio is known 
as the "quasi-neutrality" parameter and it is small when deviations from local charge 
neutrality are negligible. This is because, in most plasmas, charge imbalance can persist 
only on length scales shorter than the Debye length [13], [H] . If the latter is very small, 
the plasma is almost neutral at macroscopic scales ~ L. This requires, of course, that 
at least two types of charges, positive and negative, be present in the system under 
consideration. 

In the forthcoming sections we will show that the electric and magnetic limits can 
be obtained by making suitable assumptions on the parameters a and /3. In particular: 

• If (3 1 and a = 0(1), we obtain the electric limit; 

• If (3 <C 1 and a Cl, but a//3 = 0(1), we obtain the magnetic limit. 

3. Electric limit 

For convenience of notation, we define a smallness parameter e <C 1. Then, for the 
electric limit, we assume (5 ~ e and a = 0(1). We expands both fields in a power series 
in e, such that: E = E + eEi + . . . and B = B + eBi + . . ., where E , Ei, etc . . . are 
quantities of order unity. The density and current are also assumed to be of order unity. 
To lowest (zeroth) order in e one obtains 

V-E =p/«, (11) 

V-B = VxE = VxB = 0. (12) 
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Thus B = 2 and we can write E = — V0cb so that: 

A0 O = -p/a (13) 
Let us now proceed to first order. We find 

V • Ei = V ■ Bi = V x Ei = 0, (14) 

VxB^fl + ^V (15, 



x a dt 

which imply that Ei = 0. Introducing the vector potential at first order, Bi = V x Ai, 
and substituting into Eq. (115]) . we obtain 

V(V • Ai) - AAi = J/a - d t (V<f>o), (16) 

which reduces to 

AAi = -J/a (17) 
if we adopt the Lorentz gauge relation 

^ + V.A 1 = 0. (18) 

Putting together the results at zeroth and first order, Maxwell's equations in the 
electric limit can be written in terms of the potentials 

A0 O = - p/a, (19) 
AAi = - J/«, (20) 
or in terms of the fields 

V-E =p/a, (21) 
V • Bi = V x E = 0, (22) 

VxB 1 = i + ^. (23) 
a at 

Note that only E and B x appear in the above equations, because Ei = B = 0. 
Equations (f2Tl) - (l23l) are identical to the equations postulated by LBLL for the electric 
limit pp. As anticipated in the introduction, they can be heuristically obtained from the 
full Maxwell's equations by dropping Faraday's induction term. Here the same result 
was derived from a systematic expansion in powers of a small parameter. The present 
method also allowed us to recognize that the electric and magnetic fields actually appear 
at different orders in e. 



f Strictly speaking, the fact that V • B = V x B = implies AB = 0, but if B vanishes at 
infinity, then it must be zero everywhere. Similar considerations apply to vector fields with zero curl 
and divergence encountered later in this paper. 
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Lorentz force.— Let us now evaluate the Lorentz force per unit volume: 5F = 
pE + J x B. In our units, and using as a reference value for the force 5F = pE, 
we obtain 

5F = pE + pj x B, (24) 

where it appears that the magnetic term is of higher order. Since in the electric limit 
only E and Bi are no n- vanishing and ~ e, we get 

5F = pE + e 2 J x Bi. (25) 

Thus we conclude that first order terms in Maxwell's equations will induce a second 
order magnetic correction in the particle dynamics, which can be neglected in the present 
approximation. 

Condition of validity.— Next, we would like to derive the condition of validity for the 
electric limit as established by LBLL, i.e., |E| 3> c|B|. For this purpose, we compute 
the relative strength of the electric and magnetic fields. One can write 

|E| E E + eEi + . . . E + eE 1 + . . . . . 

- — - « = = c . 26 

|B| B B + eB 1 + ... 5 + e5i + ... V ; 

For the electric limit E\ = Bq = 0, whereas Eq ~ B\ = 0(1). Therefore we obtain: 

and since e 1, we find that |E| 3> c|B|, which is the expected result. 

Finally, the continuity equation can be obtained by taking the divergence of Eq. 
fTT5|) and using Eq. ffTT|) . which yields 

^ + V-J = (28) 
4. Magnetic limit 

This is obtained by taking (3-^-0 and a — > 0, but keeping the ratio /3/at finite. In other 
words a ~ /? ~ e < 1. By performing the same expansion as in Sec. [31 we obtain at 
zeroth order 

V ■ E = V x E = V • B = 0, (29) 

V x B = J . (30) 

Note that Gauss's law also implies that p = if we assume - as we did so far - that the 
charge density is a zeroth order quantity (a more general case will be discussed shortly). 
Thus, the magnetic limit deals with systems that are locally charge neutral, a fact that 
was already acknowledged by LBLL. 

Equation (I3"0"|) represents the magneto static limit: no free charges, only currents; no 
electric fields, only magnetic fields. We also note that the current is divergence free. 
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If we pursue the expansion to first order, we obtain 

V x Bi = V • Bi = V • Ei = (31) 
VxE 1= -^. (32) 

Summarizing the results at zeroth and first order, we can write: 

V ■ B = V • Ex = (33) 

V x B = J (34) 

VxE 1= -^, (35) 

and B x = E = 0. We note that this is the approximation of Maxwell's equations used 
in magneto-hydrodynamics (MHD). Indeed, MHD is a one- fluid theory that describes 
the motion of a fluid that carries electric currents but no electric charge. The currents 
generate a self-consistent magnetic field through Eq. (1341) . which in turn acts back 
on the fluid via the Lorentz force. Equations (133]) - (135]) are almost identical to the 
equations postulated by LBLL for the magnetic limit pQ, except that those authors 
found V ■ Ei 7^ 0. We will return on this point later. 

Introducing the vector potential Bo = V x Ao and using the Coulomb gauge 
V • A = 0, yields 

AA = -J. (36) 

Thus, it appears that the Coulomb gauge is the natural choice for the magnetic limit. 
In dimensionless units, the electric field is written in terms of the potentials: 

dA 

E = -V0-e— . (37) 
At first order, this becomes 

Ei = -V0! - (38) 
which satisfies automatically Eq. 



Charge neutrality.— The charge neutrality condition (p = 0) is correct only if one 
assumes that the density is a zeroth-order quantity in e. Different equations are obtained 
if the density is a first or second order quantity. For instance, if we assume p = ep%, 
with pi = 0(1), then Gauss's law becomes: V ■ E = p\. But in this case, we would 
have an electric field at zeroth order, which is in somewhat at odds with the spirit of a 
"magnetic" limit. 

More interestingly, we consider the case p = e 2 p2, with p 2 = 0(1). The only 
difference with respect to Eqs. ( 13~1~]) -( 132]) is that Gauss's law now reads as: V • Ei = p 2 . 
In summary, the equations at zeroth and first orders become, expressed in terms of the 
fields: 



V-Ei =p2, 
V-B =0, 



(39) 
(40) 
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V x B = J, (41) 

VxE 1= -^. (42) 

The above equations are now identical to those of LBLL for the magnetic limit. It will 
appear later that the condition p = e 2 p 2 yields the correct Lorentz transformations of 
the four-current in the magnetic limit (see Sec. [6]). 

We also point out that it was already recognized in some textbooks on classical 
electromagnetism [15] that the presence of free charges (p ^ 0) in a model based on the 
magnetic limit is a second-order effect in (3. This effect is usually neglected in standard 
MHD, which assumes that p = 0. 

Using Eq. (1381) and the Coulomb gauge condition V ■ Ao = 0, we can write Eqs. 
(}39]) - (j4"2"]) in terms of the potentials: 

AA = - J (43) 
A0! = -p 2 . (44) 

Notice that the above equations are the exact analog of the equations (fl9 |) -(l2"0 ]) that 
were obtained in the electric limit, except that the orders are inverted (electric effects 
appear at zeroth order in the electric limit, but at first order in the magnetic limit). 

Continuity equation. — Pushing the expansion to second order, we find that 

VxB 2 = ^, (45) 

so that d t (V ■ Ex) = d t p 2 = 0. This is compatible with the condition V ■ J = 0. In 
summary, the continuity equation can be written as follows: 

dtP2 + V • J = 0, (46) 

where each term is zero. 

Lorentz force.— Using Eq. (1241) and remembering that, in the magnetic limit, 
Eo = Bi = and p = e 2 p 2 , we obtain: 

5F = e 3 p2Ei + eJ x B . (47) 

Thus, there is no force at zeroth order and only a magnetic force at first order. The 
electric force appears at third order and is thus uninfluential. 

Condition of validity.— Using Eq. (1261) and considering again that, in the magnetic 
limit, Eq = Bi = while B ~ E\ = 0(1), we obtain: 

[El eE x 

]B| ~ C ~E\ ~ C€ ( ' 

and since e < 1, we have that |E| ^ c|B| . This is the condition obtained by LBLL for 
the validity of the magnetic limit. 
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5. Gauge relations 

Let us begin with the general Lorentz gauge condition 

If + V.A^O. («) 
Using our dimensionless variables, this becomes 

+ V ■ A = 0. (50) 
Now, in the electric limit we have 0i = Ao = 0, so that the gauge relation becomes 
^ + V. Al =0, (51) 

i.e., in the electric limit the natural gauge is the Lorentz gauge. 

In the magnetic limit we have 0o = Ai = 0, so that the gauge relation becomes 

6 2 ^ + V-A = 0. (52) 

At first order in e, the first term in the above equation can be neglected, so that we are 
left with: V ■ A = 0, i.e., in the magnetic limit the natural gauge is the Coulomb gauge. 

Once again, these "natural" gauge conditions, well known in the literature, arise 
automatically when applying our expansion procedure. 

6. Lorentz transformations 

In this section, we apply the above expansion technique to the Lorentz transformations 
of the space-time coordinates, the electric and magnetic fields, and the sources. We will 
obtain - in a simple and systematic way - all the known results in both limits. 

6.1. Space-time 

Under a Lorentz transformation, the space-time 4-vector (ct, x) transforms as follows: 

. v(v ■ x) . , 

x ' = x - 7 vi + 7 - i I > 53 

V 

*'=7(*-^), (54) 

where v is the relative velocity between the unprimed and primed reference frames. For 
low velocities, the Lorentz factor 7 goes like 7 ~ 1 — f3 2 /2 so that at first order in e = 
we can take 7 = 1. This yields: 

x' = x - vt, (55) 

tl =t — . 56 

c z 

Now we normalize time to T, space to L, and velocity to V = L/T. We obtain 

x' = x — vt, (57) 
t' =t-e 2 v-x. (58) 
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To second order in e, the Eq. ( 15 8 j) becomes simply £' = £, thus yielding the standard 
Galilean transformations. 

Note that, since |x| ~ L, t ~ T, we have that: |x|/(c£) ~ L/(cT) ~ e. This yields 
|x| <C c£, so that the 4-vector (ct,x) is ultra-timelike. 

6.2. Electric and magnetic fields 

The Lorentz transformations for the fields are as follows: 

E' = 7 (E + v x B) + (1 - 7 ) V ^ V - E \ (59) 

v 

1 v(v • B) 

B' = 7 (B - - v x E) + (1 - 7 ) 1 - ' . (60) 

With the approximation 7 ~ 1 (valid up to second order) and using the normalized 
units defined in Sec. HJ one obtains ((3 = e): 

E' = E + evxB, (61) 
B' = B — e v x E. (62) 

Now we expand the electric and magnetic fields in the usual way: E — Eg+eE^ + . . ., 
B = B + eBi + . . .. Substituting into Eqs. (|6T |) -( l62l) and matching order by order yields: 

E' = E (63) 

B' = B (64) 

E; = Ei + v x B , (65) 

B' 1 = Bi-vx E . (66) 

In the electric limit, Ej = B = 0, therefore the Eqs. (I63l) - (l66|) reduce to 

K = E (67) 
B; = Bx - v x E . (68) 

These are the correct Galilean transformations of the fields in the electric limit pp. 
In the magnetic limit, B x = E = 0, therefore the Eqs. (I63|) - (l66|) reduce to 

B = B (69) 
E; = Ex + v x B . (70) 

These are the correct Galilean transformations of the fields in the magnetic limit p]. 

6.3. Current and density 

The current density 4-vector (cp, J) transforms as follows: 

J' = J-7pv + ( 7 -l)^A (71) 
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Taking 7 = 1 and using dimensionless variables, we get 

J' = J - V p, (73) 
p> = p _ e 2 v . J; (74 ) 

which to first order in e yield finally 

J' = J - V p, (75) 

P = P- (76) 

The above equations are the standard Galilean transformations in the electric limit. 
Indeed, since p ~ p and |J| ~ J = Vp, one gets that cp = J/e 3> j J | - The 4-current 
(cp, J) is thus an ultra-timelike vector. 

For the magnetic limit, we assume - as was done in Sec. H]- that the density is a 
second order quantity, i.e., p = e 2 p2. Then Eqs. ( 1741) become: 

J' = J - e 2 vp 2 , (77) 

P*2 = P2 - v • J. (78) 

The second order term in Eq. fl77|) can be neglected, so that we finally obtain: 

J' = J, (79) 
p' 2 = p 2 - v ■ J. (80) 

This is the standard Galilean transformation of the current and density for the magnetic 
limit. Indeed, since |J| ~ J = Vp and p ~ p~e 2 , we have that 

- — - ~ — - ~ - (81) 
cp ce z e 

and thus | J| ~ e~ x cp ^> cp. In this case, the 4-current is an ultra-spacelike vector. Note 
also that in the strictly neutral limit (p = 0), the relation |J| ^> cp is a fortiori satisfied. 

Finally, it is possible to show [5] that the reduced Maxwell equations, both in the 
electric and in the magnetic limit, are Galilei covariant according to the respective non- 
relativistic transformations of the space-time coordinates, the electromagnetic fields and 
the sources, as derived in the preceding paragraphs. 

7. Lagrangian formulation 

In the Lorentz gauge, the Lagrangian density for an electromagnetic field obeying 
Maxwell's equations can be written in the form: 

By applying the Euler-Lagrange equations to the above Lagrangian density 
d dC ^ dC dC 



Non-relativistic limits of Maxwell's equations 



12 



(where ip is either the scalar potential or one of the components of the vector potential 
and dj is the j'th component of the gradient operator V), we obtain the usual wave 
equations for the scalar and vector potentials in the Lorentz gauge 

< 84 » 

1 d 2 A 

- AA + ^^ = ^oJ. (85) 

Let us rewrite the Lagrangian density in the usual dimensionless variables (where 
£ is normalized to p<j) ). We obtain: 

where we have defined for convenience CJ = C/a. Since we are interested in 
developments up to first order, we can disregard the two terms proportional to (3 2 . 
This yields finally: 

For the electric limit, as usual, /3 ~ e 1 and a = 0(1). Expanding the potentials 
in powers of e, one gets: £ = C'o + e£'x, with 
(V0 O ) 2 (VA ) 2 p0 o 



C' 



2 2 a 



C\ = V0 O • V0i - VA — + (89) 

a a 

Applying the Euler-Lagrange equations to O , A and their derivatives, we obtain at 
zeroth order: 



= - plot, (90) 
AA = 0. (91) 

At first order, the Euler-Lagrange equations with respect to 4>i an d Ai yield the same 
equations as above. In contrast, applying the Euler-Lagrange equations to 4>q and A 
results in 

A0i =0, (92) 
AAi = - J /a. (93) 

If the functions <p\ and Ao are zero at infinity, then they must vanish everywhere, so that 
we recover the usual equations for the electric limit expressed in terms of the potentials. 

In the magnetic limit, we have ~ a ~ e <C 1. If we also take p = e 2 p2, we obtain 
from Eq. (!H7|) at zeroth and first order 

£ ;.(M_(!M + J . Ao , (94) 
£'i = V0 O • V0i - VA • VAi - p 2 0o + J Al (95) 
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Proceeding as was done above for the electric limit, yields the equations: 

AA = - J, (96) 

A0i = -P2, (97) 



together with A0o = AAi = 0, which imply that the latter potentials vanish identically. 
Once again, we have found the correct equations for the magnetic limit. 

8. Other limits 

From the present analysis, it appears that the electric and magnetic limits are the only 
nontrivial Galilean limits of the Maxwell equations. However, they do not appear to be 
perfectly symmetric. For the electric limit, no additional assumption had to be made 
on the sources, whereas for the magnetic limit the charge density must be a second- 
or higher-order quantity in e. The magnetic limit requires that the system be "quasi- 
neutral" (using plasma physics terminology) and therefore both negative and positive 
charges must be present. 

Now, the question is whether it is possible to obtain an electric limit that is also 
quasi-neutral. In plasma physics, for instance, it is common to employ approximate 
models that are purely electrostatic (Poisson's equation) and also quasi-neutral. The 
rationale behind this approximation is that a plasma can be non-neutral only on 
distances shorter than the Debye length Xp- Over longer distances, free charges are 
screened and the plasma is basically a neutral medium, although electric fields can still 
be presently. Writing Poisson's equation in normalized units: A</> = —p/a, it is clear 
that the charge density must be of order a = (Xd/L) 2 in order to generate a finite 
electric potential. 

In the present context, one could obtain such an approximate model by first going 
to the electric limit (/3 — > 0, with a finite) and then taking a — > 0. But this procedure 
is not completely satisfactory, because there is no reason why the relative smallness 
of the parameters a and /3 should not be specified from the start. In the following 
paragraphs we suggest a way to derive, in a more rigorous way, an electric limit that is 
also quasi-neutral. 

In deriving the standard electric limit, we had assumed that /3 ~ e and a ~ 1. 
Lets us now "raise the order" by one unit, i.e., /3 ~ e 2 and a ~ e, so that we still have 
/3/a ~ e. Further, in order to have a zeroth-order electric field, we must require that 
p = ep\, with pi = 0(1). Then, Maxwell's equations for the zeroth and first order fields 
read as: 

V-E = Pl , VxB 1 = J, (98) 
V-Bi = 0, VxE = 0. (99) 

§ In one of the most popular plasma physics textbooks |13j one can read the following statement: "In a 
plasma, it is usually possible to assume n e = m and V • E ^ at the same time. This is a fundamental 
trait of plasmas, one which is difficult for the novice to understand" . This corresponds to the fact that 
a very small charge density (p <~ e) can still generate a finite electric field. 
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together with Ex = B = 0. The above equations are sometimes referred to as the 
"quasi" static" limit of Maxwell's equations, because they are obtained by neglecting all 
time derivatives in the original equations. 

In terms of the scalar and vector potential and assuming the Coulomb gauge 
V • Ai = 0, these equations can be written as: 

A0 O = -pi, AAi = -J. (100) 

From Eq. fl24l) . the Lorentz force becomes: 5F = epiE + e 3 J x B 1? and the last term 
can be neglected. Therefore, we end up with a purely electric force and a quasi-neutral 
system, because p~e~«< 1. 

One can verify that the Lorentz transformations (up to first order in e) for the 
space-time variables reduce to the standard Galilei transformations. For the fields, since 
(3 = e 2 , we have: E' = E and B' = B, and for the sources: p[ = p\ and J' = J — ewpi. 

Under the above transformations, Gauss's law in Eqs. (I9~8 l) -( l99l) is invariant, but 
Ampere's law is not. However, since the resulting force is purely electric, only Gauss's 
law V • Eo = pi needs to be taken into account. We further note that the force 
5F = epiE is also invariant under the same transformations. In summary, the resulting 
model is indeed purely electrostatic (only E counts) and quasi- neutral (since p ~ e). 

9. Conclusion 

In this work, we developed a systematic yet simple approach to the non-relativistic limits 
of Maxwell's equations. When the latter are rewritten in suitable normalized variables, 
two dimensionless parameters appear. These parameters represent the typical velocity 
normalized to the speed of light (/3) and the degree of charge neutrality of the system 
under consideration (a). 

The main result of this paper was that the non-relativistic limits of Maxwell's 
equations can be recovered by letting either or both these parameters go to zero: if 
/3 — > and a = 0(1), we recover the electric limit; in contrast, letting j3 — > and 
a — > 0, but keeping a/(3 = 0(1), leads to the magnetic limit. 

These results were obtained by expanding the electric and magnetic fields in a 
power series in the relevant smallness parameter, then matching terms at zeroth and 
first order. This procedure is both rigorous and systematic, and yields all known results 
on each of the two limits. Most known properties arise naturally within the present 
approach: for instance, the Lorentz gauge is shown to be the appropriate choice in the 
electric limit, whereas the Coulomb gauge is more suited to the magnetic limit. 

In addition, this approach revealed some hitherto overlooked subtleties of the 
theory. For instance, we learnt that: (i) the magnetic and electric fields in the limit 
equations are actually quantities at different orders in e; (ii) in the electric limit, the 
Lorentz force has no correction at first order; (iii) in the magnetic limit, the charge 
density must be at least a second-order quantity (and is usually neglected in MHD 
applications) . 
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From previous work and the present analysis, it appears that no other nontrivial 
Galilean limits exist for the Maxwell equations. We could derive an electric limit that is 
also quasi-neutral, but this limit is Galilei covariant only at zeroth order, which reduces 
to Gauss's law. 

Other limits might be found if one also expands the sources in powers of j3. When 
the Dirac four-current is expanded in powers of c _1 (i.e., (3), it displays several correction 
terms beyond the standard Schrodinger expression For instance, the charge 

density contains a second-order term, arising from the so-called Darwin correction in 
the Hamiltonian [16], which reads as: P2 = 8 ^ 2c 2 ^(^^)- In general, Galilei covariance 
will be lost at second order, but the resulting limits are still worth investigating as they 
are important to understand the coupling of the Maxwell and Dirac equations in semi- 
relativistic effective field theories [12]. The present approach should be well adapted to 
this purpose. 
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